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I.  INTRODUCTION 

This  technical  report  presents  the  analytical  solution  of  the  system  of 
coupled  ordinary  differential  equations  (ODEs)  which  describes  the  time  evolu¬ 
tion  of  an  Ideal  (l.e.,  zero  resistance)  coaxial  plasma  gun  operating  in  the 
so-called  snowplow  mode  (Refs.  1  and  2)  being  driven  by  a  single  loop  RLC  cir¬ 
cuit.  The  full  equations  appear  to  be  difficult  to  solve  analytically  in 
closed  form  although  others  (Refs.  3-13)  have  obtained  approximate,  asymptotic 
and  numerical  solutions  to  several  variants  of  the  system  of  equations.  How¬ 
ever,  In  the  special  situation  where  the  effect  of  the  coaxial  gun  on  the  cir¬ 
cuit  is  negligible  (l.e.,  the  gun  Is  fully  Influenced  by  the  circuit  but  the 
circuit  Is  negligibly  Influenced  by  the  gun)— a  state-of-affairs  called  the 
weak  coupling  limit— the  system  of  coupled  ODEs  Is  readily  solvable  In  closed 
form.  This  solution  to  the  weakly  coupled  problem  Is  presented  here,  leaving 
for  a  subsequent  report  the  details  of  the  solution  to  the  fully  coupled 
equations. 

Section  II  presents  the  physical  model  of  the  system  and  Introduces  nota¬ 
tion,  and  Section  III  describes  the  equations  governing  the  system.  Section 
IV  presents  the  solution  In  the  weak  coupling  limit,  while  Section  V  discusses 
the  validity  of  the  weak  coupling  limit.  Finally,  Section  VI  discusses  rele¬ 
vant  numerical  results. 
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II.  THE  PHYSICAL  MODEL 

A  schematic  diagram  of  the  physical  system  under  consideration  Is  shown 
In  Figure  1.  The  operation  of  this  system  Is  as  follows. 

Initially  (time  ■  0-,  Fig.  la),  the  capacitor,  with  capacitance  C,  Is 
fully  charged,  with  the  charge  on  the  positive  plate  being  q(0“)  *  Q  >  0;  the 
circuit,  with  Inherent  Inductance  LQ  and  Inherent  resistance  RQ,  Is  open;  the 
current  Is  1(0“)  •  0;  and  the  region  between  the  straight  coaxial  cylinders, 
with  Inner  and  outer  radii  r^  and  rQ,  Is  full  of  unionized  gas  of  uniform  mass 
density,  p. 

Then  (Fig.  lb),  the  switch  Is  closed  and  the  potential  difference  across 
the  Inner  and  outer  electrodes  causes  the  gas  In  a  thin  washer-like  layer  next 
to  the  rear  Insulator  to  break  down,  forming  a  conducting  plasma  sheath  there 
and  closing  the  circuit,  allowing  a  current  to  flow  and  the  capacitor  to  begin 
discharge.  Call  the  first  Instant  time  zero,  when  essentially  all  the  gas  in 
the  thin  layer  Is  In  a  plasma  state  (an  Idealization),  and  take 

q(0)  «  Q  and  1(0)  »  0  -  4(0)  (1) 

Also,  denote  the  thickness  of  the  Initial  plasma  sheath  by  a,  a  >  0. 

In  the  snowplow  mode  of  operation,  the  current  through  this  initial 
plasma  sheath  rapidly  heats  It  up,  resulting  In  the  formation  of  a  shock  wave 
that  begins  to  propagate  down  the  gun.  This  shock  wave  Initially  deposits 
some  of  Its  energy  In  a  thin  layer  of  unionized  gas  directly  in  front  of  the 
Initial  plasma  sheath,  causing  this  gas  to  Ionize  and  become  heated  plasma 
also.  Simultaneously,  the  radial  current  through  the  plasma  Interacts  with 
the  azimuthal  magnetic  field  produced  by  the  coaxial  current  to  produce  an 
axial  force  on  the  now  augmented  plasma  mass,  causing  It  to  begin  to  acceler¬ 
ate  down  the  gun.  (See  Ref.  2  for  a  more  detailed  description.)  The  net 
result  of  these  two  effects  Is  that  the  total  mass  of  plasma  in  the  gun  Is 
now  greater  than  that  In  the  original  plasma  sheath  and,  furthermore,  this 
more  massive  plasma  sheath  Is  moving  down  the  gun.  This  process  Is  repeated 
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continually  (with  the  Initial  plasma  sheath  being  replaced  by  the  updated 
plasma  sheath) ,  and  at  any  time  t  >  0,  an  accreting  mass  of  plasma  moves  down 
the  gun  with  Instantaneous  mass  equal  to  that  of  the  unionized  gas  Initially 
(l.e.,  at  t  *  0)  contained  between  the  Insulator  and  the  current  position  of 
the  leading  edge  (assumed  sharp,  flat,  and  perpendicular  to  the  gun  axis  of 
symmetry)  (Fig.  lc) .  Denote  the  distance  from  the  leading  edge  of  the  plasma 
sheath  to  the  front  of  the  rear  Insulator  by  z(t).  Note  that 

z(0)  *  A  (2a) 

Also,  take  the  velocity,  i(t),  of  the  sheath  to  be  zero  Initially: 

*(0)  *  0  (2b) 

This  definition  was  chosen  for  z  because  we  want  the  mass,  m(t),  of  the  plasma 
sheath  to  be  proportional  to  z(t). 

Also  note  that  our  definitions  of  the  circuit  quantities  Imply 
•  1(t)  «  -dq/dt  (3) 


AFWL-TR-86-137 


III.  THE  GOVERNING  EQUATIONS 

We  now  derive  the  two  equations  governing  the  time  evolution  of  the 
system:  the  circuit  equation  and  the  equation  of  motion  of  the  plasma  sheath. 

1.  THE  CIRCUIT  EQUATION 

The  circuit  Is  a  single  loop  RLC  circuit  with  capacitance  C  (other  [stray] 
capacitance  Is  neglected);  with  Inductance 

L(t)  -  lQ  +  Lg(t)  (4) 

where  Lg(t)  Is  the  time-dependent  Inductance  of  the  portion  of  the  gun 
between  the  front  of  the  rear  Insulator  and  the  leading  edge  of  the  plasma 
sheath : 


Lg(t)  *  C(V2w)ln(ro/r1^*^  (MKS)  (5) 

Lnote  that  Lg(o)  *  A(y0/2ir)ln(r0/r^ )  >  0];  and  with  resistance 
R(t)  *  Rq  +  Rw(t)  +  Rp(t),  where  Rw(t)  is  the  resistance  contributed  by  the 
portion  of  the  gun  walls  lying  between  the  front  of  the  rear  Insulator  and 
the  leading  edge  of  the  plasma  sheath  (the  resistance  of  the  portion  of  the 
walls  bounding  the  rear  Insulator  is  Included  In  RQ),  and  Rp(t)  is  the 
resistance  contributed  by  the  plasma.  Assume  here  that  the  gun  is  ideal, 
l.e.,  Rw(t)  *  0,  and  that  the  plasma  has  Infinite  conductivity,  so  that 
Rp\t)  «  0;  hence,  R(t)  ■  RQ.  (Alternatively,  one  may  assume  Rp(t)  Is  some 
nonzero  constant  which  has  been  absorbed  into  RQ.) 

The  circuit  equation  is  standard  except  for  the  presence  of  a  variable 
Inductance.  Note  how  this  variable  Inductance  Is  handled  in  both  the  Kirchhoff 
Law  treatment  and  the  energy  conservation  treatment  of  the  circuit  (Fig.  Id). 


a.  Kirchhoff  Law 


We  have 

(q/c)  -  Rq1  -  L0(d1/dt)  +  (Vb-Va)  -  0 

where  Vfa-Va  Is  the  Induced  EW ,  5,  across  Lg.  But 

Vva  *  5  *  -d*(t)/dt  -  -d[Lg(t)1]/dt  »  -Lg(t)(d1/dt)  -  1(dLg/dt)  (6) 

where  *(t)  Is  the  flux  linking  the  variable  inductor,  and  the  inductance  Lg(t) 
is  defined  by 

♦(t)  *  Lg(t)1(t)  (7) 


(q/c)  -  [R0  +  (dLg/dt)]1  -  [Lq  +  Lg(t)] (di/dt)  -  0 
Now  recalling  that  i  ■  -q  (Eq.  3), 

D-0  +  +  [R0  +  *Lg(t)M  +  (1/C)q  »  0  (8) 

•  • 

which  is  the  desired  circuit  equation.  Note  that  when  Lg  >0,  Lg  is 
equivalent  to  a  resistance,  and  the  11  term  is  in  some  sense  a  dissipative 
one;  this  is  made  more  precise  in  the  energy  conservation  treatment. 

b.  Energy  Conservation 

To  use  this  approach,  one  needs  expressions  for  the  amount  of  energy 
stored  in  the  capacitor  and  in  the  inductor  at  any  time  t.  For  an  ideal  capac¬ 
itor  with  constant  capacitance  C,  this  is  considered  to  be  the  amount  of  energy 
invested  in  bringing  the  capacitor  to  charge  q(t)  from  charge  zero,  namely, 
(l/2C)q2(t);  and  for  an  ideal  inductor  with  constant  inductance  lQ,  this  is 
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considered  to  be  the  amount  of  energy  Invested  in  bringing  the  inductor  to 
current  i(t)  from  current  zero,  namely  (l/2)L0i2(t) . 

For  an  ideal  variable  inductor,  the  energy  stored  is  also  considered 
to  be  the  energy  invested  and  is  calculated  as  follows:  Let  W(t)  denote  the 
amount  of  energy  invested  in  "currenting"  the  variable  inductor,  with  induc¬ 
tance  Lg(t),  from  time  zero  to  time  t  and  let  c(t)  denote  the  induced  EMF 
across  this  inductor  at  time  t.  Then 


According  to  Eqs.  9  and  10,  one  may  interpret  (l/2)L(t) i 2  as  the  amount 
of  energy  stored  in  the  inductor  at  time  t  and  further  available  as  electro¬ 
magnetic  energy  to  drive  the  circuit,  and  (1/2) i  (t')L(t')dt‘  as  the  amount 
of  energy  stored  in  the  inductor  at  time  t  and  not  further  available  as 
electromagnetic  energy  to  drive  the  circuit,  with  (1/2) i 2L  being  the  rate  of 
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“dissipation"  of  ircuit  electromagnetic  energy  by  the  variable  inductor. 

Of  course,  this  Is  not  dissipation  in  the  usual  sense  because  the  lost 
electromagnetic  energy  all  appears  as  (mechanical)  kinetic  energy  of  the 
deforming  variable  inductor;  in  our  case,  It  appears  as  kinetic  energy  of  the 
plasma  sheath. 

2.  EQUATION  OF  MOTION  OF  THE  PLASMA  SHEATH 

The  motion  of  the  plasma  sheath  is  that  of  an  accreting  mass,  with  instan¬ 
taneous  momentum  m(t)z(t),  being  pushed  down  the  gun  by  the  force  associated 
with  the  energy  "dissipated"  in  the  variable  inductorT  By  our  definition  of 
z(t),  the  mass  is  given  by 


where  p  is  the  density  of  the  unionized  gas.  The  model  makes  no  precise 
statement  about  the  thickness  of  the  plasma  sheath  and,  hence,  none  about  its 
density. 

An  expression  for  the  force  on  the  plasma  sheath  may  be  derived  using 
either  of  the  two  points  of  view  following: 

a.  All  the  energy  dissipated  in  the  variable  inductor  is  transformed 
into  translational  kinetic  energy  of  the  sheath;  that  is,  the  sheath  is  con¬ 
sidered  to  be  mechanically  rigid.  (For  this  reason,  the  sheath  is  often  called 
a  slug.)  If  F  denotes  the  force  on  the  sheath,  then  ( 1/2) i 2Lg  *  F*z(t).  Now 
by  Eq.  5, 


*Lg  *  C(M0/2ir)ln(r0/ri)]z(t) 
so  whenever  z(t)  t  0 

F  *  C(p0/4ir)ln(ro/ri)]q2  (11) 


a 
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b.  All  of  the  energy  dissipated  In  the  variable  Inductor  Is  expended  In 
the  work  done  by  the  magnetic  field  pressure  In  pushing  the  sheath  down  the 
gun;  that  Is,  there  is  no  field  diffusion  into  the  plasma.  This  Is  consistent 
with  the  fact  that  we  have  taken  Rp  *  0.  The  pressure  on  the  rear  of  the 
sheath  at  radius  r  Is  B2(r)/2yQ,  where  B(r)  *  yQ1/2irr,  so  that 


F  ■  (2ir/2y0)  /  B  (r)rdr  «  [(y^iOln^/r^]  q* 


in  agreement  with  Eq.  11, 


The  equation  of  motion  Is  then 

P7r(r02  "  rl)at  {zi)  ’  (V4ir)ln(ro/r1}  * 


ZZ  +  Z 


i2  *  [(u0/4^p)1n(r0/ri)/  (rQ2  -  ri2)]^ 


In  summary,  the  governing  equations  for  our  physical  system  are 
(Eqs.  1,  2a,  2b,  5,  8  and  12): 


[Lq  +  bz(t)]q  +  [RQ  +  bz(t)]q  +  (1/C)q  =  0 


2  .2 


••  «  •  «  . 
zz  +  z  =  aq  (t) 


where 


q (0 )  «  Q  >  0  ${0)  *  0  z(0)  =  A  >  0  z(0)  =  0 


2  2 


a  *  (yQ/4ir  p)ln(r0/ri)/(r0  -  ri  ) 


b  *  (u0/2n)ln(r0/ri) 
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We  also  require  that 

z(t)  >  0  for  all  t  >  0  (18) 

The  parameters  a  and  b  may  be  thought  of  as  coupling  parameters:  A  Is  a 
measure  of  the  Influence  of  the  circuit  on  the  gun  and  b  Is  a  measure  of  the 
Influence  of  the  gun  on  the  circuit.  Our  weak  coupling  limit  Is  (approximately) 
the  one  In  which  the  parameter  b  Is  considered  to  be  small.  This  should  not 
be  confused  with  the  small  coupling  limit  (Ref.  8)  In  which  the  parameter  a  is 
considered  to  be  small. 
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IV.  SOLUTION  IN  THE  WEAK  COUPLING  LIMIT 


w: 


The  details  of  the  solution  of  the  governing  equations  In  the  weak  coupling 
(WC)  limit  follow.  In  the  next  section,  the  conditions  are  examined  under 
which  the  WC  limit  Is  applicable. 

Given  $(t),  Eq.  14  Is  easily  solved  for  arbitrary  Initial  conditions  by 
noting  that  It  may  be  written  as 

2  2  2  2 
d  (z  )/dt  *  aq 


Then 


:(t)  ■  ±[z(0)  +  2z(0)z(0)t  +  2a  fo  Jo  '’!(u)dud<]i 


(19) 


For  our  conditions  we  have 


:(t)  »  a[l  +  (2a/&  )ff  q  (u)dudv]* 

*/o  *'o 


(20) 


Of  course,  q  depends  upon  z  via  Eq.  13  so  q(t)  Is  not  In  fact  "given"  in 
general.  However,  in  the  special  case  that 


bz(t)  «  L 


(21a) 


and 


bz(t)  «  Rq  (21b) 

for  all  times,  t,  of  interest,  the  circuit  equation  becomes 

L0q  +  Roq  +  (1/C)q  -  0  (22) 

which  equation  Is  easily  solvable  for  q(t);  Insertion  of  this  q(t)  into 
Eq.  20  then  yields  the  time  evolution,  z(t),  of  the  plasma  sheath.  Equations 
21  and  22  give  the  precise  definition  of  the  WC  limit. 

11 
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Concerning  the  solution  for  z(t),  note  fro*  Eq.  20  that 


x(t) 


.2 

q  (u)du 


(23) 


so  that 

z(t)  >  0  for  all  t  >  0  (24) 

Also,  note  fro*  Eq.  20  that 

z(t)  £  A  >  0  for  all  t  >  0  (25) 

Hence,  the  plasma  sheath  Is  always  moving  down  the  gun  (and  never  moving 
backwards  up  the  gun).  Since  the  gun  has  finite  length,  the  plasma  sheath 
eventually  reaches  the  muzzle  end,  say  at  time  ■  tg.  The  times  of  Interest 
referred  to  after  Eqs.  21  are  then  0  <  t  <  tg.  The  model  has  no  applicability 
for  t  >  tg. 

The  solution  to  the  WC  limit  circuit  equation  (Eq.  22),  with  Initial 
conditions  on  q  and  q  given  by  Eq.  15,  Is  well  known  to  be.  In  the  non- 
crltlcally  damped  case, 

q(t)  •  Q  exp(-RQt/2l0)[cos  fit  +  (R0/2nLQ)s1n  fit]  (26) 


where 


a  .  t(i/i0c)  -  0>oV«lo')]*  *  o 

and  where.  If  Q  *  1«  with  u  real,  then  cos  fit  *  cos(iwt)  «  cosh  wt  and 
sin  fit  ■  sln(lwt)  ■  Islnh  wt.  In  the  critically  damped  case, 

q(t)  -  Q  exp(-R0t/2L0)[l  ♦  (R0/2LQ)t]  (27) 

Of  course,  the  general  Initial  value  problem  with  q(0)  arbitrary  could  be 
solved  easily,  but  for  notatlonal  simplicity  later  on  [aside  from  the  Issue 

12 
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of  th«  natural  Initial  conditions  for  our  physical  nodal ],  we  choose  not  to  do 
so:  tha  solution  for  z(t)  with  4(0)  •  0  will  already  be  notatlonally  cumber¬ 
some.  The  sane  comment  applies  to  tha  Initial  value  problem  for  z(t):  *e 
use  £q.  20  rather  than  Eq.  19  to  find  z(t)  from  q(t).  Substituting  Eq.  26  or 
27  Into  Eq.  20,  one  gets,  respectively  (using  v  ■  Rp/2l0  and  6  ■  RQ/2DL0), 

-l(t)  •  41  ♦  (aO'/EA*)  ^e“2fint^(l  -  4*)cos  2Qt  *  4  sin  2ot 

♦  £fs*  4  (i/4* )]  4  i  j  4  Ja  4  (i/4)jot  -  ija  4  (i/al)A|  (28) 


12  2  /  -2vt  2  t  ^ 

1  4  (aQ  /2a  )  le  [(vt)  ♦  2wt  4  ij)  4  vt  -  ^  j  (29) 

These  are  the  desired  solutions  for  the  notion  of  the  plasma  sheath  in  the 
WC  limit. 

For  completeness,  we  also  present  expressions  for  the  velocity,  z(t), 
in  the  two  cases.  They  are 

z (t)  -  (nA)(A/2)(aQa/4A,)(1^4*)|e"2aQt[4cos  2flt-sin  Zflt-[4*(l/4)]Jt(;/6)  (30) 


2  I  I  -2v  t  2  I 

z(t)  •  (vA)(A/z)(aQ  /4a  )|l-e  [2(vt)  4  2vt  4  1] j  (31) 

Note  that  the  Issue  of  the  qualitative  dependence  of  the*-solut1on  of  the 
circuit  equation  (Eq.  13)  upon  the  coefficients  of  that  equation  has  been 
Ignored.  That  Is,  for  any  fixed  time  t  the  nature  of  the  solution  to  Eq.  13 
depends  upon  the  sign  of 


»rf*ere 


D(t)  •  [R0  -  (4l0/C)]  4  bA(t) 


A(t)  •  bilt)  4  2R0i(t)  -  (4/C)z(t) 


•’.V/AV'-T-*’ 
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Using  Eqs.  13-15  It  my  be  shown  that  A(0)  »  -4a/C,  A(0)  >  0  ■  A(0),  and 
A(0)  *  4aR0Q*/AC*(L0  ♦  bA)a  >  0  so  that  A(t)  Is  not  constant  in  tlm.  One 
ray  therefore  envision  a  choice  of  RQ,  LQ,  C  such  that  C(t)  changes  sign 
at  least  once  during  [o,tg],  and  so  speculate  that  the  nature  of  the  damping 
changes  during  the  motion  of  the  sheath. 


1 

1 

♦ 

1 
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V.  VALIDITY  OF  THE  WEAK  COUPLING  LIMIT 

In  this  section  some  conditions  are  examined  under  which  the  WC  approxl 
nation  Is  applicable.  The  length  of  the  coaxial  gun  gas  compartment  Is 
denoted  by  Zg. 

First,  consider  the  requirement  Imposed  by  Eq.  21a.  An  upper  bound  on 
bz(t)  for  0  £  t  <  tg  Is  bZg.  If  we  take  x  «  y  to  mean  ax  *  y  (usually 
a  2.  10),  then  It  Is  necessary  and  sufficient  that 

Lo  i  obZg  *  «(w0/2*)1n(r0/ri>Zg  (32) 

Next,  consider  the  requirement  Imposed  by  Eq.  21b.  An  upper  bound  on 
bi(t)  for  0  £  t  <  tg  Is  not  as  readily  obtainable  as  that  on  bz(t);  the 
result  Is 

bi(t)  x  ba*q/(L(JC)i.  0<t<tg  (33) 

which  we  will  demonstrate  shortly.  It  Is  then  sufficient  that 

Rq  i  aba*Q/(L0C)*  (34) 

If  VQ  Is  the  Initial  potential  difference  across  the  capacitor,  then 
Q  ■  CVQ  and  so  Eq.  34  becomes 

«0  *  •b**Vc/L0)* 

Combining  this  with  Eq.  32,  the  condition  on  RQ  then  becomes 

*  <“»VC/V* 

•  (»#*c)  ln(r0/r1)/^Z9(r01  -  r,1^  (3S) 

Noting  that  uQ  *  4n  x  10" 7  H/m  (MKS) ,  we  get  sufficient  conditions  on  LQ 
and  Rq,  In  terms  of  the  other  system  parameters,  for  the  WC  limit  to  be 
valid: 

15 
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L0  >  0.2aZgln(ro/ri)  (tiH) 


R0  i  7.979  x  10-s  a*(v02C^  ln^/r^/j^Z^2  (mo)  (37) 

We  now  prove  the  result  In  Eq.  33.  First*  It  follows  from  Eq.  10 
(energy  conservation)  and  the  condition  L  >  0  that,  for  0  <  t  <  tg,  we  have 

^012(t)  <  |t(t)12(t)  *  (1/2C)Q2 

so  that 

q2(t)  s.  Q2/LqC,  Oititg  (38) 

Hence,  frjm  Eqs.  23*25,  for  0  <  t  <  tg,  we  have 

z(t)  i  [»/z(tn(Q,/L0c)tJ  <  (./i)^VL0cjts 

Since  z  Is  bounded  above  on  compact  [o,tg].  It  attains  a  global  maximum 

there.  This  maximum  Is  attained  at  an  Interior  stationary  point,  trttere 

z  »  0,  or  at  the  endpoint  t  »  t  ,  where  z  >.  0  [for  If  z(tg)  <  0  then  z(tg)  > 

z(tg)  which  contradicts  the  fact  that  z(tg)  Is  a  global  maximum].  In  any 

case.  If  we  denote  by  t  a  time  at  which  z  achieves  Its  global  maximum, 

max 

and  denote  z(tMJ  by  i  .  then  z( tMV)  >  0  and  Eqs.  14  and  18  Imply  that 
max  max  max 


4  ^  *[<(w] 


Now  using  Eq.  38  In  Eq.  39  gives 


z2  <  aQ2/L  C 
max  -  v  '  o 


whereupon  Eq.  33  follows  lamed lately. 


iByiyijvv 
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VI.  NUMERICAL  RESULTS 

In  this  final  section,  numbers  and  graphs  relevant  to  the  preceding 
theory  are  presented.  We  arbitrarily  choose  a  *  10;  upon  absorption  of  a*, 
the  constant  In  Eq.  37  becomes  2.523  x  10“*. 

Tables  1  and  2  give  the  minimum  values  of  LQ  required  by  Eq.  36  and  of 
Rq  required  by  Eq.  37,  respectively,  for  various  representative  values  of 
r0/ri  and  Zg.  In  the  latter  case,  VQ,  r^  and  C/p  are  fixed  at  the  typical 
values  5  x  10*  V,  0.03  m  and  10"a  F/kg-m*  respectively;  this  last  choice  is 
consistent  with  the  range  of  values  of  C  and  p  representative  of  plasma  guns 
operating  In  the  snowplow  mode,  namely,  10“*  <  C  <  10“*  (F)  and  10"5  £  p  <. 
10'*  (kg/m*). 

Figures  2-4  show  1(t),  z(t)  and  i(t)  for  representative  underdamped, 
critically  damped  and  overdamped  cases,  respectively.  For  these  graphs  we 
use  the  dimensionless  forms 

1/oQ  -  (1  +  «l)  e"acts1n  nt 

from  Eq.  26,  z/ A  from  Eq.  28  and  i /PA  from  Eq.  30  for  the  noncrltlcal ly  damped 
case,  and  we  use  the  dimensionless  forms 

i/\>Q  »  vte’vt 

from  Eq.  27,  z/A  from  Eq.  29  and  z/va  from  Eq.  31  for  the  critically  damped 
case. 


Concerning  the  choices  of  values  of  S  and  dimensionless  coupling 
strength  (see  Eqs.  28-31) 

<  =  aQ2/A* 

note  the  following: 
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TABLE  1.  MINIMUM  VALUES  OF  L0  (uH)  FOR  WHICH 
THE  WC  LIMIT  IS  VALID  FOR  VARIOUS 
VALUES  OF  Z  (m)  AND  r/r1 ;  a  -  10 


x?/r< 

VN 

1.1 

2 

10 

0.1 

0.019 

0.139 

0.461 

0.5 

0.096 

0.695 

2.31 

1 

0.191 

1.39 

4.61 

TABLE  2.  MINIMUM  VALUES  OF  Rq  (ma)a  FOR  WHICH 
THE  WC  LIMIT  IS  VALID  FOR  VARIOUS 
VALUES  OF  Zg  (m)  AND  rQ/r^;  C/p  «  10"J 
F/kg-ffls,  Vp  *  5  x  10'  V,  r1  »  0.03  m; 
a  ■  10 


V»,r( 

V\ 

1.1 

2 

10 

0.1 

27.7 

53.2 

30.8 

0.5 

12.3 

23.7 

13.7 

l 

8.73 

16.8 

9.73 

aValues  of  L  given  in  Table  1,  namely  L^n,  are  assumed  here.  For 


values,  Lq,  of  inductance  greater  than  lgin,  the  minimum  values  of  R 
given  here  decrease  as  (Lg1n/L0)*  {see  Eq.  34). 
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a.  The  choices  6  *  0.0119  for  the  underdamped  case  and  |«[  *  1.407 

for  the  overdamped  case  (where  by  our  convention,  6  is  purely  imaginary) 
correspond  to  LQ  *  335  nH,  C  *  3  yF  and  RQ  =  8  mft  or  RQ  =  950  «Q,  respectively 
(with  n  »  1.003  x  10s s"1  or  u  *  1.008  x  106s"1,  respectively).  (For  critical 

damping  with  this  LQ  and  C,  RQ  =  668j  mfi  and  v  =  0.997  x  106s  l.)  Tables  1 

and  2  show  that  these  choices  of  L0  and  RQ  are  sufficient  for  the  validity 

of  the  WC  limit  under  some  circumstances,  e.g.,  when  rQ/r^  *  1.1  and  Zg  < 

2  m  (and  ri  =  0.03  m  and  VQ  *  5  x  10“  V).  These  choices  of  RQ,  L0  and  C  also 

make  n,  u>,  v  *  1.0  x  10*s_1  so  that  the  numerical  values  of  at,  «t  and  vt  on 

plot  abscissas  are  very  nearly  equal  to  the  numerical  values  of  t,  in  units 
of  microseconds,  to  which  they  correspond. 

b.  The  possible  sensible  values  of  x  span  a  range  of  at  least  twelve 
orders  of  magnitude,  from  1.83  x  10“7/A2  to  1.45  x  10S/A2,  corresponding  to 
the  two  sets  of  values  {rQ/r^  *  10,  r-j  =  0*05  m,  C  =  1  yF,  p  *  10~3kg/m3, 

V  =  2.5  x  103  V)  and  {r Jr.  =  1.1,  r.  =  0.01  m,  C  =  1000  yF,  p  *  10-5kg/m3, 

VQ  =  10s  V},  respectively.  The  actual  value  of  A  is  unimportant  since,  as 
can  be  seen  from  Eqs.  28-31,  it  just  establishes  a  basic  scale  length  in  terms 
of  which  the  quantities  z  and  i  are  measured.  For  definiteness  we  take  A  = 
10’3  m.  Then  1  x  10*1  <  ic  <  1  x  1010.  Clearly,  for  small  x,  z(t)  is  rela- 

'v 

tively  small  and  z(t)  is  relatively  flat.  To  exhibit  some  structure  in  the 
graphs  we  choose  a  large  value  of  x,  namely  x  =  1  x  103 ;  this  corresponds 
to  the  set  of  values  rQ/r^  *  1.1,  r..  *  0.03  m,  C  =  3  yF,  p  =  3  x  10  kg/m3, 
and  VQ  =  5  x  10“  V. 
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